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Abstract. We calculate the asymptotic merit factor, under all cyclic rotations 
of rows and columns, of two families of binary two-dimensional arrays derived 
from the quadratic character. The arrays in these families have size p X q, 
where p and q are not necessarily distinct odd primes, and can be considered 
as two-dimensional generalisations of a Legendre sequence. The asymptotic 
values of the merit factor of the two families are generally different, although 
the maximum asymptotic merit factor, taken over all cyclic rotations of rows 
and columns, equals 36/13 for both families. These are the first non-trivial 
theoretical results for the asymptotic merit factor of families of truly two- 
dimensional binary arrays. 



1. Introduction 

We consider an array of size n x m to be an infinite matrix A — {uij) of real- 
valued elements satisfying 

ttij = unless < * < and < j < m. 

The array is called binary if Uij takes values only in {1,-1} for all i,j satisfying 
< i < n and < j < m, and is called ternary if aij takes values only in {0, 1, — 1}. 
Given integers u and v, the aperiodic autocorrelation of ^ = (oij) at displacement 
(m, v) is defined to be 

Ca(u,w) := aijOj+uj+v 

We refer to an array of size n x 1 as a sequence of length n, abbreviating the array 
(flio) to (ai) and the aperiodic autocorrelation Ca{u,0) to Cyi(u). 

Binary arrays with small out-of-phase aperiodic autocorrelation have a wide 
range of applications in digital communications and storage systems, including 
radar [1] and steganography [23]. Ideally, we would like to find a binary array 
A of size n X m satisfying 

(1) \CAiu,v)\<l for all (0,0), 

in which case, A is called a Barker array [Tj. However, it was recently shown by 
Davis, Jedwab, and Smith [5] that a (truly two-dimensional) Barker array must 
have size 2x2. (Barker sequences of length n, namely n x 1 Barker arrays, are 
known for n S {2, 3, 4, 5, 7, 11, 13}, and any other Barker sequence must have even 
length [22 greater than lO^^ [TSj.) 
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Since the Barker array criterion ([T]) is too restrictive for array dimensions ex- 
ceeding 2 X 2, it is natural to define a measure for the collective smallness of the 
aperiodic autocorrelation values of a binary array. One such measure is the merit 
factor, which is defined for a binary array A = (uij) of size n x m with nm > 1 to 
be 

F(A) ■= (^"^)' 

Let Fn^m denote the maximum value of F{A) taken over all 2"™ binary arrays A 
of size nx m, and abbreviate i^„,i to F„. We note that the mean of 1/F{A), taken 
over all binary sequences A of length n, equals 1 — 1/n [19]. The argument of [19] 
easily generalises to two dimensions: the mean of 1/F{A), taken over all 2"™ binary 
arrays A of size n x m, equals 1 — l/(nm). It follows that Fn.m > nm/{nm — 1), 
which asymptotically equals 1 and provides a first benchmark result. 

A number of theoretical and computational results on Fn are known (see [llj for a 
survey). One line of research is to calculate F„ for small values of n. At present, F„ 
has been calculated for all n < 60 (see |TT, Fig. 1], for example). The largest values 
of Fn currently known are F13 = ^ ~ 14.08 and Fn = which are attained 
by Barker sequences. The computational analysis of Fn quickly becomes infeasible 
as n grows. Another line of research is therefore to construct particular infinite 
families of binary sequences of increasing length and to calculate their asymptotic 
merit factor. 

The only non-trivial infinite families of binary sequences for which the asymp- 
totic value of the merit factor is known are Rudin-Shapiro sequences [T7] , Legendre 
sequences [3], and m-sequences [Uj, together with some generalisations of these 
families [10] , [15] , [20] , [12] , [13] . The largest proven asymptotic merit factor of a bi- 
nary sequence family is 6, which is attained by cyclically rotated Legendre sequences 



(see Theorem 2.1 1. There is also considerable numerical evidence, though currently 
no proof, that an asymptotic merit factor greater than 6.34 can be achieved for a 
family of binary sequences related to Legendre sequences [31 . 

Much less is known about the value of Fn^m for n,m > 1. Eggers [T] Tab. 4.2] 
computed Fn^m for nm < 21 and found lower bounds on Fn^m for nm < 121. 
Although the data supplied in [7] are very limited, it is apparent that F„ „i tends 
to be smaller than Fnm- The largest value of Fn^m for n,m > 1 reported in [7] 
equals F4 4 = ^ ~ 5.33. However, an elementary construction technique gives 
binary arrays with larger merit factor. Given two sequences A = (a^) and B = {bj) 
of length n and m, respectively, we follow [S] in defining the product array Ax B := 
(aibj). A straightforward calculation shows that 

(2) CAxBiu,v)^CAiu)CB{v), 

from which we deduce 

F{AxB) = {^^F{A))i^^F{B] 



- 1. 



Let A and B be Barker sequences of length 13 and 11, respectively. It follows 
from ([3| that F{A x A) ~ 6.80, F{A x B) ~ 6.27, and F{B x S) ~ 5.81. Another 
consequence of ^ is 

F F 
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Currently, no theoretical results on the asymptotic merit factor of families of bi- 
nary truly two-dimensional arrays are known. Bomer and Antweiler [5] analysed the 
merit factor of several binary array families numerically. Among the investigated 
families, two types of array families related to the quadratic character appeared to 
have largest merit factor. The arrays in the first family were proposed by Calabro 
and Wolf [5| and have size p y. q, and the arrays in the second family were pro- 
posed by Bomer, Antweiler, and Schotten [3_ and have size p x p, where p and q 
are (not necessarily distinct) odd primes. Both families can be considered as two- 
dimensional generalisations of Legendre sequences. The authors of [2] successively 
applied three operations, namely rotations of rows and columns, stairlike rotations 
of rows and columns, and proper decimations, and computed the maximum value 
of the merit factor for arrays of small sizes taken from these families. They then 
remarked [2, p. 8] that 

". . . , for large arrays, the ACF [aperiodic autocorrelation function] 
merit factors of both classes [the above mentioned array families of 
square size] appear to tend to 3.", 

and asked for a theoretical explanation of this observation. 

In this paper we study the merit factor of two families of binary arrays. The 
arrays in the first family, called Legendre arrays, have size p x q, where p and q 
are (not necessarily distinct) odd primes, and contain as a special case the arrays 
proposed by Calabro and Wolf [5 . The arrays in the second family, called quadratic- 
residue arrays, have size p x p, where p is an odd prime, and contain as a special 
case the arrays proposed by Bomer, Antweiler, and Schotten [3 . We calculate, 
under certain conditions on the growth rate of p relative to q, the asymptotic merit 
factor at all rotations of rows and columns for both array families. In particular, 
we show that for both families the asymptotic merit factor equals || ~ 2.77 for an 
optimal rotation of rows and columns. Although we only maximise the merit factor 
with respect to the first operation considered in [2], namely rotations of rows and 
columns, this result does not support the conclusion of |3] quoted above. For all 
other (non-optimal) rotations of rows and columns, the asymptotic merit factor of 
quadratic-residue arrays is larger than that of square Legendre arrays. 



Given an odd prime p and a positive integer m, let GF(p™) be the finite field 
containing elements. Whenever convenient, we treat integers after reduction 
modulo p as elements in GF(p). The quadratic character o/GF(p™) is the function 
X : GF(p™) ^ M defined by 



2. Two Families of Binary Arrays 








+1 



1 



for a = 



for a not a square in GF(p'") 
otherwise. 



This function is multiplicative: 



(4) 



X(a)x(^) = X{ab). 
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If m = 1, then {a\p) :— x(a) is the Legendre symbol. A Legendre sequence L = [li) 
of prime length p > 2 is defined by 

^ [l for i = 
' [ (* |p) for 1 < i < p. 

If the initial element in a Legendre sequence is changed to zero, so that 

= {i\p) for < i < p, 

then we call L a ternary Legendre sequence. 

In what follows, we present two families of binary arrays, which can be considered 
as two-dimensional generalisations of Legendre sequences. Let p and q be two (not 
necessarily distinct) odd primes, and let Vp,g be the set of ternary arrays V — {vij) 
of size p X q satisfying 

f 1 for (i = and < j < g) or {j = and < i < p) 
I otherwise. 

Then Vp^q contains 2^'+'^"^ arrays, each having p + q — 1 nonzero elements. Given 
ternary Legendre sequences L and K of length p and q, respectively, we define 
a Legendre array X of size p x g to be a binary array of size p x q that can be 
written as 

X = L X K + V for some V G Vp,g. 
For example, the array X — (xij) of size p x q, given by 

-1 for j = and < i < p 

+1 for 1 = and 1 < j < q 

(* \ p)U \q) ioT 1 < i < p and 1 < j < <?, 

is a Legendre array of size p x q. This particular array was originally defined by 
Calabro and Wolf [5] , and its merit factor properties were investigated numerically 
in P] . In the original paper [5] such an array was called a "quadratic-residue array" . 
We use the term Legendre array to distinguish it from our second family of binary 
arrays. 

Let p be an odd prime, let x be the quadratic character of GF(p^), and let {a, a'} 
be a basis for GF(p^) over GF(p). Following Bomer, Antweiler, and Schotten [5, 
we define a quadratic-residue array Y = (jjij) of size p x p to be a binary array of 
size p X p satisfying 



+1 or —1 for i = j = 

X{ia + ja') for 0<i,j <p, ^ (0,0). 



The class of quadratic-residue arrays {yij) satisfying yoo = +1 was defined by 
Bomer, Antweiler, and Schotten [S], and its merit factor properties were investigated 
numerically in [2] . In our analysis it will be convenient to change the leading element 
in a quadratic-residue array to zero. Accordingly, we define the ternary quadratic- 
residue array of size p x p to be the ternary array Z = (zij) of size p x p given 

by 

(5) := x(ia + ja') for 0<i,j <p. 
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Next we define an operation acting on an array to produce a new array of the 
same size. Given an array A = (a^j) of size n x m and real numbers s and t, the 
rotation Ag^t is the array B — (bij) of size n x m given by 

(6) &ij = a(i+LnsJ)modn, 0+L™tJ)modm for < i < n and < j < TO. 

If A is a sequence of length n, we abbreviate Ag^o to As- 

The asymptotic merit factor of a Legendre sequence was calculated for all rota- 
tions by H0holdt and Jensen |!9j. 



Theorem 2.1 (H0holdt and Jensen 9J). Let L be the Legendre sequence of prime 

1 



length p > 2, and let s be a real number satisfying — ^ < s < ^. Then 



hm FiL^^^'^'^-^' 

p— >-oo 



The constraint — | < s < | in Theorem 



2.1 



is for notational convenience only 
since by definition Ag is the same as A^^i for every sequence A and all real s. The 
maximum asymptotic merit factor of a rotated Legendre sequence Lg is 6, which 
occurs for s = ±i. 



3. Calculation of the Merit Factor of an Array 
Given a positive integer n, let 

be a primitive (2n)th root of unity. Let A — (aij) be an array of size n x m. The 
generating function of A is defined to be the power series 

If A is a sequence of length n, we write A{x) for A{x, y). 

The next lemma shows how the merit factor of A can be computed from the val- 
ues A(Cni Cm)- This approach generalises to two dimensions the method of Il0holdt 
and Jensen [S] to compute the asymptotic merit factor of a sequence. 

Lemma 3.1. Let A be an array of size n x m. Then 

2n-l 2m-l 



u.v i—0 j—0 

Proof. Straightforward manipulation shows that 

A{x, y)A{x'\y'^) = ^ Ca(u, for a; ^ and y ^ 0, 

U,V 

and therefore, 

(7) \Aix,y)\' =J2CA{u,v)x-"y-- for |a;| = |y| = L 
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Using this identity, an elementary calculation gives 

2n-12m-l 2ri-12r?i-l 



^ I O J. ^ I I L _l_ ^ ^ I L J. ^ t t I, _l_ 

E l^(c.a)r = ^E E iMQXDWicu^' 

u,v 

as required. □ 

4. The Merit Factor of Legendre Arrays 

In this section we compute the asymptotic merit factor of Legendre arrays for all 
rotations, subject to certain conditions on the growth rate of the dimensions. We 
first record a result on the aperiodic autocorrelation of rotated ternary Legendre 
sequences, which arises as an immediate corollary of |201 Thm. 3]. 

Proposition 4.1. Let L be the ternary Legendre sequence of prime length p > 2, 
and let s be a real number satisfying — ^ < s < |. Then, as p oo, 

4 E [^^^^ HI' = i + 8 (1*1 - if + oip-'i^ogpf)- 



We note that, as explained after [20, Thm. 3], we can recover Theorem 2.1 



from Proposition |4.1[ We also need the following bound for the magnitude of a 
polynomial over C at a (2ci)th root of unity, in terms of its magnitudes at dth roots 
of unity. 

Lemma 4.2. Let d > 1 be odd, and let Ad C[x] have degree at most d — 1. Then 
\A{Ci)\ < (21ogd) max |^(Cf )| for integer j. 

Q<k<d 

Proof. The bound is trivial in the case that j is even. We may therefore take j to 
be odd, writing j — 2£ + d for some integer £ so that = — Cj^. It is then sufficient 
to bound Now by La grange interpolation we have 

^(-) = ^E|^cf A(cf), 

k=0 ^d- 



and so, since d is odd. 



1 2 

|^(-Cr)| < J E \/-2e , /■2fe| l^(C 

" k=0 l^rf ~^^d I 



f) 



2 '^^^ 1 

< - max Ufd*")! -75— r. 

- d o<k<d I ^^'^ ^ 1 + Cw™ 

The result follows from the inequality 'Y^m-tii^lX^ + Cd"N — ^logd (which holds 
since d is odd, see [HI p. 625], for example). □ 

The next theorem gives, under certain conditions on the growth rate of p relative 
to g, the asymptotic merit factor of all 2^*+*"^ Legendre arrays of size p x g for all 
rotations. 
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Theorem 4.3. Let he an infinite set of products of two not necessarily distinct 
odd primes, and let N take values only in M. Write N = pq for odd primes p 
and q, and suppose that 



(8) 



q p 

— and — ^ as N ^ oo. 
p2 qZ 



Let X be a Legendre array of size p x q, and let s and t be real numbers satisfying 
— I < s,t < ^. Then 



(9) 



lim F{Xs,t) 



l + 8{\s\^iy l + 8{\t\-iy 



- 1. 



N 



Proof. Let L and K be the ternary Legendre sequences of length p and q, re- 
spectively, and write T -.^ L x K. Notice that Tg^t = Lg x Kt and Ts^t{x,y) = 
i^(x)i4:t(2/). Thenfrom ([2]) 

U.V U V 

The condition Q implies that p and q grow without bound as ^ oo. We can 
therefore apply Proposition 4.1 to give 
(10) 

Y.^CTUu,v)f^\l+8{\s\-\f] [l+8(|t|-i)'l +o(l) asTV^^. 



1 



{pq)' 

Define 

(11) 

We claim that 
(12) 



A(7V) 



[pqf 



A(iV)^0 asiV^cx). 



The theorem then follows from ( [lO| and ( 11 ) and the fact Cx^ t (0, 0) — pq. 

It remains to prove the claim p2|). By the definition of a Legendre array, there 
exists an array V G ^p,q such that X^^t — Ts,t + ^s,t- From Lemma 3.1 we then 
have 



A(iV) = 



1 



4(pg)s 



2p-l 2g-l 



2p-l 2q-l 



E E |7^M(c;,q) + v;,t(c;,q)r- E E I^M(c;,q)| 

i=0 j=0 1=0 j=0 



^ 2p-12g-l 

(13) < E E 



Now for a, 6 e C the identity 

\a + 6|4 = |a|4 + + 4[3fJ(a6)]' + 2 |ap • l&p + 4 jap • 3fi(a6) + 4 |6|2 • 3fJ(a&) 
gives the inequality 

\a + - |a|4 < 4 |a|3 • |6| + 6 |ap • |6p + 4 |a| • + 
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Apply this bound to ( 13 1 to obtain 

2p-12g-l 



+ 



(14) 



3 
1 
1 

4(w)' 



1=0 j=0 
2p-l 2g-l 

E E I^M(C,q)r-|^..t(c,q)| 

1=0 j=0 
2p-l 2q-l 

i=0 j=Q 
2p-l 2q-l 

rE EI^^.*(C'C^')r- 



i=0 j=0 



Given a ternary Legendre sequence A of length d, it is well known (see [HJ p. 182], 
for example) that |y4(Cj'^)| < d}/'^ for each integer k. It is easily verified that 
this implies |j4r(Cj'^)| < d}^"^ for each integer k and all real r. Therefore, since 
Ts.t{x,y) = Ls{x)Kt{y), Lemma 



4.2 



gives 



\TsMi^ Cg) I < ^{pqf^ log(p + q) for aU integers i and j. 
Substitute into ( 14 1 to give 

(15) A(Ar) < 256 + ^^^' ^1+96 ^^^^^^^ '^^^' s. + 16 ^3- 



(w)i/2 



(pg)3/2 



r'S'4, 



where 



2p-12ij-l 

^^-i^E EI^M(C^a)r- 

■'^^ i=0 j=0 

From ([7]) and a straightforward calculation we obtain 

52 = Cy^,(0,0) 

= p + (7 - 1. 

The Cauchy-Schwarz inequality gives Si < {82)^^^, and since \Vs,t{x, y)\ < p+q — l 
for \x\ = |t/| = 1, we also have Si < {p + q— 1)^^^S'2 for £ > 2. Substitution into ( 15 ) 
gives 

1/2 



A{N) < 256 
+ 16 



pq 

(p + g-l)2 



(log(p + q)y 



96 



p + q-l 



pq 



(logip + q)y 



log(p + q) 



(p + g - 1) 

(pg)2 



(pg)3/2 

Now, using the condition (Isl), we readily verify our claim (12). 



□ 



4.3 



from the restriction — i < s,t < i 



There is no loss of generality in Theorem 
since As^t is the same as As+i.t and As^t+i for every array A and all real s and i 
We note that the condition ([8| can be relaxed for particular Legendre arrays. For 
example, let L and K be the Legendre sequences of length p and q, respectively. 
Then X = L x K is a, Legendre array. From ^ and Theorem 2.1 we conclude 
that Q holds under the relaxed condition p — > cx) and g — > 00 as — > 00. 
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5. The Merit Factor of Quadratic-Residue Arrays 

In this section our goal is to calculate the asymptotic merit factor of a quadratic- 
residue array of size px p at all rotations. We shall assume throughout this section 
that p is an odd prime. Write the pth roots of unity as 

ej e^/^^'":i/p for integer j. 

Then, since p is odd, we have 

{Cp-0<i<2p}^ {e, : < j < p} U {-e^ : < j < p}. 

Therefore, given an array A of size p x p, Lemma [3T] asserts that 

(16) j2^CA{u,v)r 

= ^2 E (l^(e.,e,)l' + l^(-Q,e,)l' + l^(e.,-e,)|^ + |^(-Q,-e,)n. 

0<i,j<p 

Our objective is to find an asymptotic expression for the sum on the right-hand side 



of the identity (16), where A is a rotated ternary quadratic-residue array. Since 
a ternary quadratic-residue array and a quadratic-residue array differ in only one 
element, this will be sufficient to compute the asymptotic merit factor of a rotated 



quadratic-residue array. Before we analyse the sum in (16), we shall need several 
technical results, which we state in the next subsection. 

5.1. Auxiliary Results. The following lemma evaluates the generating function 
of a ternary quadratic-residue array at pth roots of unity. 

Lemma 5.1. Let x be the quadratic character o/GF(p^), and let Z be a ternary 
quadratic-residue array of size p x p, as defined in ([s]) . Then there exists a basis 
{/3,/3'} for GF(p2) over GF{p) such that 

Z{ek,ei>) = {-l)'^px{k(3 + ep') for all integers kj. 
Proof. Let Tr : GF(p^) ^ GF{p) be the trace function given by 

Tr(a;) = a; -|- x^, 

and for b e GF(p^), let ipb ■ GF(p^) ^ C be the additive character of GF(p^) given 

by 

Mx) :=e^2.Tr(M/f>. 

It is readily verified that 

(17) Tr(ax-|-6?/) = aTr(a;) + 6Tr(j/) for a, & € GF(p). 

We choose {/3, 13'} such that {a, a'} (appearing in the definition ([s]) of Z) and 
{(3, (3'} are dual bases, that is, 

(18) Tr(a/?) = 1, Tr(a/3') = 0, Tr(a'/3) ^ 0, Tr(a'/3') = 1. 

Such a basis is guaranteed to exist |161 p. 58]. Given integers i, j, fc, £, we then have 



by (17) and (18) 



Tii{ia+ja'){k(3 + ei3')) = ik + ji, 
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and therefore 

0<2,J<p 

aGGF(p2) 

by putting a :— ia + ja' . The above sum is cahed a Gaussian sum, and it is well 
known that 

^ x{a)Ma)^{-lY-^PX{b) 

aeGF(p2) 

(see [ini PP- 199-201], for example). This proves the lemma. □ 
Our next lemma bounds a certain character sum and evaluates it in special cases. 
Lemma 5.2. Let x be the quadratic character of GF{p^), and define 
(19) r!(K,A,M):= J2 

x{x)x{x + k)x{x + X)x{x + for k, X, fi e GF{p'^) 

xeGF(p^) 

and 
(20) 

/ := {(k, k, 0) : k e GF{p^)} U {(k, 0, k) : k G GF{p^)} U {(0, k, k) : k e GF{p^)}. 
Then 

p^ — 1 for K = 
p2 _ 2 for K^O 



(21 



and 
(22) 



pin, X, fi)\ < 3p for (k, X, ^i) ^ I . 
Proof. Since x is multiplicative by Q, 

fi(«,K,0)= ixi^ix + K))]^ 

xeGF{p^) 
2;eGF(p2)\{o,-K} 

using x(0) = and [x(a;)]^ = 1 for all nonzero x £ GF(p^). This proves (21). 

To prove (22 1, we use a special case of a result on multiplicative character sums 
with polynomial arguments [16, Thm. 5.41], which can be stated as follows. If 
/ £ GF(p^)[a;] is a monic polynomial of positive degree d that is not a square (that 
is, f{x) cannot be written as f{x) = [g{x)]'^ for some polynomial g £ GF(p^)[a;]), 
then 



(23) 



< {d-l)p. 



xeGF(p^) 

For all (k, A, fi) ^ /, the polynomial f{x) :— x{x + k)(x + X){x + /i) is not a square. 
Using the multiplicativity (|4| of x, application of (231 gives (22). □ 

The following technical lemma can be obtained from the results of [9^ and [TS] . 
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Lemma 5.3. Define 

p-i ^ 

(24) r(fc, ^, to) := > r for integer k,£,m. 

^ (1 + e»)(efe + e^)(e^ + e.)(em + e.) 



Then 



jor k = V (mod p) 



(25) r(fc,/c,0) 



48 

p2 1 



— for fc ^ (mod p) 



au- 



di 



(26) \nkJ,m)\<{p\ogpY 

Vl<k,l,m<p 



Proof. The identity ( 25 1 was established in [9l p. 162, Cases 4 and 5] . The bound ( [26] 
follows from the inequahty YTi=o + ^jI — plogP (see [El p. 625], for example) 
since 



|r(fc,^,TO)|< 

0<k,e,m<p 0<i,k,i,m<p 



el 



/p-i 



(1 + ei)(efe + ei){ee + e^){em + e^) 



^11 + 

\i=0 ' ■/ 

5.2. Asymptotic Merit Factor Calculation. We are now in a position to anal- 
yse asymptotic behaviour of the sum on the right-hand side of the identity (16), 
where A is a rotated ternary quadratic-residue array. We split the analysis into the 
following three lemmas. 

Lemma 5.4. Let Z be a ternary quadratic-residue array of size p x p, and let s 
and t be real numbers satisfying — ^ < s,t < ^. Then, as p ~> oo, 

^ Y \Z.,M,e,)\'^\+0{p-'). 

0<i,j<P 

Proof. Lemma |5 . 1 1 implies 



I for i = j = (mod p) 
I p otherwise. 



Then, using the easily verified identity 

(27) Z,,,(e.,e,) = 67Lf^J,-L^'*Jz(e„e,), 

we find that 

0<i,J<P ^ ^ 

as required. □ 
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Lemma 5.5. Let Z he a ternary quadratic-residue array of size p x p, and let s 
and t be real numbers satisfying — ^ < s,t < |. Then, as p oo, 



J2 |^s,t(-e„6,)r = i+4(|s|-i)2 + 0(p-i(logp)4), 



0<ij<p 



^ J2 KM,-^j)\'^i+m-\r+o{p-\\ogpr). 

Proof. Given an array A, let A'^ denote the transpose of A. Since (Zs^t)'^ = [Z'^^t^s 
and Z^ is again a ternary quadratic-residue array, it is sufficient to prove the first 
statement in the lemma. 

Let P g C[x] be a polynomial of degree p — 1, and let i and j be integer. We 
shall make use of the Lagrange interpolation formula 

(28) 



n ^ — ^ 



k=0 



Efc + 



(see [SJ Eq. (2.5)], for example). It follows that 

2 Efe 

*(-£», Cj) = - ^Zs,t{<^k,£j) 

^ k=0 



ffc + 



Set S := [psj and T := [pij, and use (27) and Lemma 5.1 to obtain 



Z,,t(-e„e,) - 2(-l)Ve-^^e-^x(fc/? + j/?' 



fc=0 



Cfc + 



where x is the quadratic character of GF(p^) and {/?,/?'} is some basis for GF(p^) 
over GF(p). Since we also have 

p-i 

Zs,ti-e,, £,) = 2(-l)4ieJ 'kXikP + jP') 



k=0 



e/c + e?: ' 



we find that 



|Z,,t(-e„e,)|^ = 16 E ef-a4-cX(a^+j/3')x(&/3+J/3')x(c/3+J^')xW+j^') 

0<a,fc,c,d<p 

ea ev 6f^ £7" 



Use the definition (24) of r{k,£,m) to write 
p-i 

(29) V — - - = 

^ Eq + £(, + €i Ec + Ed + 



fa + £6 + ei ec + + 



Ec-a r(6 - a,c - a,d - a), 



so that after variable relabeling 



^ |Z^,f(-ei,ej)| =16 ef-^+mQ 

-i— 0<a,k,£,m<p 

X{ap + i/3')x(a/3 + J/?' + mx{ap + jp' + ^/3)x(a/3 + iP' + m/3) r(fc, I, m). 
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Put X := a/3 + jf3' and note that af3 + j/3' ranges over GF(p^) as a and j range from 
to p — 1. By the definition ( 19 1 of ^{k, A, /i) we therefore obtain 

^ |Z,,t(-e„e,)f-16 4-i+.n^t^[kfi,tp,mfi)T{k,t,m), 
Let the set / be as defined in ([20l), and write 



(30) 
where 



4p' 



0<ij<p 



B = 



4 



^ikl3, iP. m/3) r(fc, i, m) 

0<k,£,m<p 



0<fc,£,m<p 



Using (22 1, the sum A can be bounded as 

12 



0</c,£,m<p 



< 12p-^(logp)^ 



by (26). Therefore, 



(31) ^ = 0(p-^(logp)'') asp-^oo. 

Using symmetry of fl{K,X,^) and T{k,£,m) with respect to interchanging their 
arguments, the sum B can be written as 

p-i 



(32) 



B = ^ r!(0, 0, 0) r(0, O, O) + ^ J2 ^(^^^ 0)r(A:, fc, O) (26,. + e 



2S\ 
k ) 



AppUcation of (21 ) and ( |25| ) to evaluate D,{kf3, kl3, 0) and r(A:, fc, 0) then gives 

^'24 
II 



(p2 + 2)(p2-i) 2(p2-2)^'2 + ef-i 



12p4 

We wish to apply the identity 



k=l 



(33) Er 



2 2 



24 



for integer j satisfying |j| < p 



(see [HI p. 621], for example). Since ^ < s < i, we have 125" — 1\ < p for all 
sufficiently large p, which allows us to apply ( 33 ) to ( 32 ) to obtain, for all sufficiently 
large p, 



B 



(p2 + 2)(j52-l) (p2_2)2 p2_2/|25-l| 1 



12p4 



4p4 



By the definition of S, we have S ^ ps + 0(1) as p — !■ oo, and therefore. 



(34) 



B =\+A{\s\-\f + 0[p-^) asp^oo. 



The proof is completed by substituting (31) and (34) in (30). 



□ 
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Lemma 5.6. Let Z be a ternary quadratic-residue array of size p x p, and let s 
and t be real numbers satisfying — ^ < s,t < ^. Then, as p oo, 



4p' 



^ J2 \Z.A~^.,~e,)\'^l+64{\s\-iri\t\-\f + 0{p-\logpf 



0<i.j<p 



Proof. The idea of the proof is similar to that of the proof of Lemma |5.5| The 
main difference is that we now have to apply interpolation of Zg_t{x,y) in both 
indeterminates. 

Let i and j be integer. By the interpolation formula (28) we have 



P 



k=0 



Efe + 



Apply the interpolation formula ( 28 1 again to obtain 
4 



P 



o<k,e<p 



Efc + Ei Q + Cj 



Set S := [ps\ and T [pt\. Then by (27) and Lemma 5.1 we get 



P o<M<p + Q + 

where x is the quadratic character of GF(p^) and {/?,/?'} is some basis for GF(p^) 
over GF{p). Since we also have 

we obtain 



E E 



, , , 

6— a+d— c b' —a -\-d —c 



0<a.b,c.d<p Q<a' ^b' ,c.' ,d' -Cp 

X{afi + a'/3')x(6/3 + b' fi')x{cfi + fi')x[dfi + 



Ea + Ej Ca' + £6 + Ef,' + Cc + Cc' + Q + Cd' + 



Use (29), relabel the summation indices, and use the definition (19) of Q{k,X,ii) to 
give 



0<i j'<p 



E E 



S T 
^k-l+m '^k'-l'+m' 



0<fe,£,m<p 0<fe',£',m'<p 

r2(fc/3 + k'P',l(3 + f /?', m/3 + m'l3') T{k, £, m) T{k',£', m'). 



Let / be the set defined in (|20l, and write 



(35) 



Ap' 



^ E \Zs,t{-^^,-ej)\' = A + B, 



0<i,j<p 
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where 



Q<k,e,r)i,k' ,e' ,7n' <p 
(kP + k' fi' ,eil+£' ,ml3 + m' II' )ei 



n{kf3 + k'f3', ip + ep', m/3 + m'/3') r(fc, I, m) r(fc', f , m') 

64 



Q<k,£,m,k' ,i' ,m' <p 



n{kp + (.p + f + m'P') T{k, m) r(fc', f , m'). 



From ( 22 1 we have the bound 

192 



^ \ 0</c/, m<p / 



A\ < 

< 192p"^(logp) 



by (126|), giving 

(36) A = 0{p-^{logpf) asp^oo. 

We use symmetry of n{K, X, fi) and T{k,£,m) with respect to interchanging their 
arguments to partition the sum B further as 

(37) B = Bi + B2 + B3 + ^4, 
where 

64 

Si = ^r!(o,o,o) [r(o,o,o)]2 

64 

^2 = ^ E 0) 0) r(0, 0, 0)(2efe + ef ) 

^ fe=i 

64 

^3 - ^ E ^'Z^'' 0) r(0, 0, 0) Tik', k', 0){2ek' + e^J) 

^ k' = l 

64 

Bi = —, ^ikf3 + k'(3',kl3 + k'l3',0)T{k,k,0)rik',k',0)i2ekek'+efeir). 

^ l<k,k'<p 



Apphcation of (21) and (25) to evaluate fl{K,K,0) and r{k,k,0) then gives 

^ (p2 + 2)2(p2_i) 
36 p6 



2 p4 _ 4 2 + e 



3 |l-e,P 

o „4 4 o -I- 
3 3 p6 A. |1_,,|2 

P^-2 2 + 



U Z; V A A 

B4 = 16^—^ > — 

„6 1 
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We complete the proof by evaluating the asymptotic behavior of the terms Bi, B2, 

-B3, and 

The term Bi: The term Bi becomes 



(38) 



5i = ^ + 0(p- 



as p — )■ 00. 



The terms B2 and B^: Since — | < s < ^ implies that 125* — 1| < p for all 
sufficiently large p, we can use the identity ( 33 ) to write 

2' 



B2 



2 p4 



p 



2 



|2^-1| 
P 



Then, since S — ps + 0{l) as p co, we obtain 



(39) B2 = ^ + ^{\s\- + 0{p-') asp ^00. 
We proceed similarly for the term B3 and find that 

(40) B3 = ^ + l{\t\-\y+0{p-') asp ^00. 
The term B4: The term B4 can be rewritten as 



B, 



16 



p- 



p" 



/p-i 



1 



/p-l 25-1 



,2T-1 



fife I 



efcl 



16^ 



Noting that 125* — 1| < p and \2T — 1| < p for all sufficiently large p, 
application of ( 33 ) gives 

2 (p2-2)(p2-l)2 

9 p6 
-2 2 



p° 



p2 

2 



\2S-1\ 
P 



24 



|2r- 11 



24 



(41) 



Since S = ps + 0(1) and T — pt + 0{1) as p ^ 00, we finally obtain 
B, = l+ [8i\s\ - i)2 - i] [8{\t\ - If - i] + O(p-i) as p ^ 00. 



The result now follows by substituting the asymptotic forms ( 38 ) , ( 39 1 , ( 40 ) , and ( 41 1 
into ([37| and then (ISTI) and ((Sel) into ((SSl). 



□ 



We are now able to prove the main result of this section. 



Theorem 5.7. Let Y be a quadratic-residue array of size px p, and let s and t be 

real numbers satisfying —h<s,t<-k. Then 



lim F{Y,,t) 

p—^oo 



''{\t\ 



Proof. Let Z he a ternary quadratic-residue array of size p x p. From ( 16 ) we have 



0<i,j<p 



0<ij<p 



4p' 



0<'id<p 



0<i,j<p 
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Using Lemmas |5.4[ |5.5[ and |5.6| to determine the asymptotic behavior of the sums 
on the right-hand side, we obtain 



(42) \J2i(^z.A^,v)r 

u.v 



''{\t\ 



O {p ^(logp)®) as p ^ oo. 



Now, since the binary array Y differs from the ternary array Z only in a single 
position, wc have 

Cy^ = Cz.^u.v) + S{u,v), 

where 



(43) 



\Siu,v)\ < 



1 for —p < u,v < p 



otherwise. 
Then, by the Cauchy-Schwarz inequality, 



1/2 / ^ 1/2 

2 \ ( V^r;;/„, „,m2 



so that by (43) and (42) 



1 



Y,[CY.A^,v)r~J2[Cz.An,v)]' 



= 0{p ^) as p — ?> oo. 
The theorem follows from (|42]) by noting that Cy.^tiO, 0) = p^. □ 



6. Concluding Remarks 

We have computed the asymptotic value of the merit factor of Legendre arrays 
(under certain conditions on the growth rate of their dimensions) and of quadratic- 
residue arrays, for all rotations of rows and columns. The asymptotic merit factor 
of rotated Legendre arrays and rotated quadratic-residue arrays is shown in Fig- 
ures [l] and [2j respectively. The maximum asymptotic merit factor, taken over all 
rotations, equals H ~ 2.77 for both array families, which occurs at the rotations 
(s,i), where s,t & {|, |}. However, at all other rotations, the asymptotic merit 
factor of quadratic-residue arrays is larger than that of Legendre arrays. On the 
other hand, an advantage of Legendre arrays is that they are not restricted to be 
square. 

In |4j , the authors exhibited a family of sequences A obtained by appending an 
initial fraction of a rotated Legendre sequence to itself. Based on partial expla- 
nations and extensive numerical computations, it was conjectured in U that this 
sequence family has asymptotic merit factor greater than 6.34. Under the assump- 
tion that this conjecture is correct, the corresponding square product array A x A 
has asymptotic merit factor greater than 2.93, by (|3|. This suggests that the max- 
imum asymptotic value of the merit factor of the two array families considered in 
this paper, namely || ~ 2.77, can be surpassed by another array family. 
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Figure 1. The asymptotic value of F{Xs^t) (under conditions 
specified in Theorem 4.3), where X is a Legendre array. 




Figure 2. The asymptotic value of F{Ys^t), where F is a 
quadratic-residue array. 



In closing, we remark that instead of studying the merit factor of two-dimensional 
arrays of size nxm, Gulliver and Parker |8J studied the (suitably generalised) merit 
factor of d-dimensional arrays of size 2x2x---x2. In [8], the merit factor of several 
families of such arrays was calculated. In particular, the largest asymptotic value, 
as d — cx), of the merit factor of a family of d-dimensional arrays considered in [S] 
equals 3. 
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